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Abstract 

We consider a network composed of two interfering point-to-point links where the two transmitters can exploit 
one common relay node to improve their individual transmission rate. Communications are assumed to be multi- 
band and transmitters are assumed to selfishly allocate their resources to optimize their individual transmission rate. 
The main objective of this paper is to show that this conflicting situation (modeled by a non-cooperative game) has 
some stable outcomes, namely Nash equilibria. This result is proved for three different types of relaying protocols: 
decode-and-forward, estimate-and-forward, and amplify-and-forward. We provide additional results on the problems 
of uniqueness, efficiency of the equilibrium, and convergence of a best-response based dynamics to the equilibrium. 
These issues are analyzed in a special case of the amplify-and-forward protocol and illustrated by simulations in 
general. 

Index Terms 

Cognitive radio, game theory, interference channel, interference relay channel, Nash equilibrium, power allocation 
games, relay channel. 



I. Introduction 

A possible way to improve the performance in terms of range, transmission rate or quality of a network 
composed of mutual interfering independent source-destination links, is to add some relaying nodes in the 
network. This approach can be relevant in both wired and wireless networks. For example, it can be desirable 
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and even necessary to improve the performance of the (wired) link between the digital subscriber line (DSL) 
access multiplexors (or central office) and customers' facilities and/or the (wireless) links between some 
access points and their respective receivers (personal computers, laptops, etc). The mentioned scenarios give 
a strong motivation for studying the following system composed of two transmitters communicating with 
their respective receivers and which can use a relay node. The channel model used to analyze this type of 
network has been called the interference relay channel (IRC) in flU where the authors introduce a channel 
with two transmitters, two receivers, and one relay, all of them operating in the same frequency band. The 
main contribution of HAH is to derive achievable transmission rate regions for Gaussian IRCs assuming 
that the relay is implementing the decode-and-forward protocol (DF) and dirty paper coding. 

In this paper, we consider multi-band interference relay channels and three different types of protocols at 
the relay, namely DF, estimate-and-forward (EF), and amplify-and-forward (AF). One of our main objectives 
is to study the corresponding power allocation (PA) problems at the transmitters. To this end, we proceed in 
two main steps. First, we provide achievable transmission rates for single-band Gaussian IRCs when DF, EF, 
and AF are respectively assumed. Second, we use these results to analyze the properties of the transmission 
rates for the multi-band case. In the multi-band case, we assume that the transmitters are decision makers 
that can freely choose their own resource allocation policies while selfishly maximizing their transmission 
rates. This resource allocation problem can be modeled as a static non-cooperative game. The closest works 
concerning the game-theoretic approach we adopt here seem to be (HO Elm and ||9llfT0llfTTi . In J6]|, the 
authors study the frequency selective and the parallel interference channels and provide sufficient conditions 
on the channel gains that ensure the existence and uniqueness of the Nash equilibrium (NE) and convergence 
of iterative water-filling algorithms. These conditions have been further refined in 0. In J9]|, a traffic game 
in parallel relay networks is considered where each source chooses its power allocation policy to minimize 
a certain cost function. The price of anarchy lfT2l is analyzed in such a scenario. In [10], a quite similar 
analysis is conducted for multi-hop networks. In ifTTl . the authors consider a special case of the Gaussian 
IRC where there are no direct links between the sources and destinations and there are two dedicated relays 
(one for each source-destination pair) implementing DF. The power allocation game consists in sharing the 
user's power between the source and relay transmission. The existence, uniqueness of, and convergence to a 
NE issues are addressed. In the present paper however, we mainly focus on the existence issue of an NE in 
the games under study, which is already a non-trivial problem. The uniqueness, efficiency, and the design of 
convergent distributed power allocation algorithms are studied only in a special case and the generalization 
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is left as very useful extension of the present paper. 

This paper is structured as follows. Sec. HI describes the system model and assumptions in multi-band 
IRCs. Sec. Un] provides achievable transmission rates for single-band IRCs. These rates are exploited further 
in multi-band IRCs (as users' utility functions) analyzed in Sec. ITVl where the existence issue of NE in the 
non-cooperative power allocation game is studied. Three relaying protocols are considered: DF, EF, and AF. 
Sec. |IV] provides additional results on uniqueness of NE and convergence to NE for the AF protocol. Sec. 
IVl illustrates simulations highlighting the importance of optimally locating the relay and the efficiency of 
the possible NE. We conclude with summarizing remarks and possible extensions in Sec. |VI] 

II. System Model 

The system under investigation is represented in Fig. [IJ It is composed of two source nodes Si, S2 (also 
called transmitters), transmitting their private messages to their respective destination nodes T>i, V 2 (also 
called receivers). To this end, each source can exploit Q non-overlapping frequency bands (the notation 
(q) will be used to refer to band q G {1, . . . , Q}) which are assumed to be of unit bandwidth. The signals 
transmitted by Si and S 2 in band (q), denoted by x[ 9 ' and X^\ respectively, are assumed to be independent 
and power constrained: 

Q 

We {1,2}, J2 E \ x i q) \ 2 < p i- (D 

9=1 

For i G {1,2}, we denote by O^f 1 the fraction of power that is used by Si for transmitting in band (q) that 
is, E|xf' ? ' ) | 2 = 9^ Pi. Additionally, we assume that there exists a multi-band relay 1Z. With these notations, 
the signals received by T>i, V 2 , and 1Z in band (q) express as: 
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where Z\ ~ M(Q,N- ), i G {l,2,r}, represents the Gaussian complex noise on band (q) and, for all 
G {1, 2} 2 , h\f is the channel gain between «Sj and T>j and h\l is the channel gain between <Sj and 1Z in 
band (q). The channel gains are considered to be static. In wireless networks, this would amount, for instance, 
to considering a realistic situation where only large scale propagation effects can be taken into account by 
the transmitters to optimize their rates. The proposed approach can be applied to other types of channel 
models. Concerning channel state information (CSI), we will always assume coherent communications 
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for each transmitter-receiver pair (Si,T>i) whereas, at the transmitters, the information assumptions will be 
context dependent. The single-user decoding (SUD) will always be assumed at T>i and V 2 . This is a realistic 
assumption in a framework where devices communicate in an a priori uncoordinated manner. At the relay, 
the implemented reception scheme will depend on the protocol assumed. The expressions of the signals 
transmitted by the relay, Xr q \ q e {1, Q}, depend on the relay protocol assumed and will therefore also 
be explained in the corresponding sections. So far, we have not mentioned any power constraint on the 
signals Xr . Note that the signal model (O is sufficiently general for addressing two important scenarios. 
If one imposes an overall power constraint Y^q=i ^<\Xr 9 ^\ 2 < P r , multi-carrier IRCs with a single relay can 
be studied. On the other hand, if one imposes ElX^] 2 < Pr Q \q G {1, ■■■,Q}, multi-band IRCs where a 
relay is available on each band (the relays are not necessarily co-located) can be studied. In this paper, for 
simplicity reasons and as a first step towards solving the general problem (where both source and relaying 
nodes optimize their PA policies) we will assume that the relay implements a fixed power allocation policy 
between the Q available bands (E|#| 2 = P^ q \ q e {1, Q}). 

To conclude this section, we will mention and justify one additional assumption. As in ll4lll3l lfT3l . the 
relay will be assumed to operate in the full-duplex mode. Mathematically, it is known from [14] that the 
achievability proofs for the full-duplex case can be almost directly applied the half-duplex case. But this 
is not our main motivation. Our main motivation is that, in some communication scenarios, the full-duplex 
assumption is realistic (see e.g., [Q]| where the transmit and receive radio-frequency parts are not co-located) 
and even more suited. In the scenario of DSL systems mentioned in Sec. H the relay is connected to 
the source and destination through wired links. This allows the implementation of full-duplex repeaters, 
amplifiers, or digital relays. The same comment can be applied to optical communications. 

Notational conventions 

The capacity function for complex signals is denoted by C(x) = log 2 (l + x); for all a E [0,1], the 
quantity a stands for a = 1 — a; the notation —i means that — % = 1 if i = 2 and — % — 2 if % = 1; for all 
complex numbers c G C, c*, |c|, IZe(c) and Im(c) denote the complex conjugate, modulus and the real and 
imaginary parts respectively. 

III. Achievable transmission rates for single-band IRCs 

This section provides preliminary results regarding the achievable rate regions for the IRCs assuming DF, 
EF, and AF protocols. They are necessary to express transmission rates in the multi-band case. Thus, we do 
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not aim at improving available rate regions for IRCs as in |fT3l and related works lQ31l[T6lllfT7i In the latter 
references, the authors consider some special cases of the discrete IRC and derive rate regions based on the 
DF protocol and different coding-decoding schemes. In what follows, we make some suboptimal choices for 
the used coding-decoding schemes and relaying protocols which are motivated by a decentralized framework 
where each destination does not know the codebook used by the other destination. This approach, facilitates 
the deployment of relays since the receivers do not need to be modified. In particular, this explains why we 
do not exploit techniques like rate- splitting or successive interference cancellation. As we assume single- 
band IRCs, we have that Q — 1. For the sake of clarity, we omit the superscript (1) from the different 
quantities used e.g., X^ becomes in this section X;. 

A. Transmission rates for the DF protocol 

One of the purposes of this section is to state a corollary from 0. Indeed, the given result corresponds 
to the special case of the rate region derived in [3 J where each source sends to its respective destination a 
private message only (and not both public and private messages as in 0). The reason for providing this 
region here is threefold: it is necessary for the multi-band case, it is used in the simulation part to establish 
a comparison between the different relaying protocols under consideration in this paper, and it makes the 
paper sufficiently self-contained. The principle of the DF protocol is detailed in 031 and we give here 
only the main idea behind it. Consider a Gaussian relay channel where the source-relay link has a better 
quality than the source-destination link. From each message intended for the destination, the source builds 
a coarse and a fine message. With these two messages, the source superposes two codewords. The rates 
associated with these codewords (or messages) are such that the relay can reliably decode both of them 
while the destination can only decode the coarse message. After decoding this message, the destination 
can subtract the corresponding signal and try to decode the fine message. To help the destination to do so, 
the relay cooperates with the source by sending some information about the fine message. Mathematically, 
this translates as follows. The signal transmitted by Si is structured as Xj = X i0 + y^-pX^. The signals 
Xj and X ri are independent and correspond to the coarse and fine messages respectively; the parameter 
h>i represents the fraction of transmit power the relay allocates to user i, hence we have V\ + u 2 < 1; the 
parameter Tj represents the fraction of transmit power <Sj allocates to the cooperation signal (conveying the 
fine message). Therefore, we have the following result. 

Corollary 3.1 (S^U): When DF is assumed, the region in (TJ]) is achievable; for i £ {1, 2}, where j = —i, 
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\h^nPi \ c ( \hg\ 2 Pi + \h ri \ 2 iSjP r + 2TZe(h ii h* ri )^r l P l u l P r \ \ 
h jr \ 2 rjPj + N r )' \\h ji \ 2 P j + \h r tfvjP r + 2Ke{h ji h* n )- s /T j P j v j P r + Nij] 

(z/i, zaj) G [0, l] 2 5.?. i/j + ^ 2 < 1 and (n, r 2 ) G [0, l] 2 , n+r 2 < 1. 

In a context of decentralized networks, each source Si has to optimize the parameter in order to 
maximize its transmission rate In the rate region above, one can observe that this choice is not independent 
of the choice of the other source. Therefore, each source finds its optimal strategy by optimizing its rate 
w.r.t. t*(tj). In order to do that, each source has to make some assumptions on the value tj used by the 
other source. This is precisely a non-cooperative game where each player makes some assumptions on the 
other player's behavior and maximizes its own utility. Interestingly, we see that, even in the single-band case, 
the DF protocol introduces a power allocation game through the parameter representing the cooperation 
degree between the source Si and relay. In this paper, for obvious reasons of space, we will restrict our 
attention to the case where the cooperation degrees are fixed. In other words, in the multi-band scenario, 
the transmitter strategy will consist in choosing only the power allocation policy over the available bands. 
For more details on the game induced by the cooperation degrees the reader is referred to 0. 

B. Transmission rates for the EF protocol 

Here, we consider a second main class of relaying protocols, namely the estimate-and-forward protocol. A 
well-known property of the EF protocol for the relay channel lfl4l is that it always improves the performance 
of the receiver w.r.t. the case without relay (in contrast with DF protocols which can degrade the performance 
of the point-to-point link). The principle of the EF protocol for the standard relay channel is that the relay 
sends an approximated version of its observation signal to the receiver. More precisely, from an information- 
theoretic point of view |[T4|. the relay compresses its observation in the Wyner-Ziv manner |[T8l . i.e., knowing 
that the destination also receives a direct signal from the source that is correlated with the signal to be 
compressed. The compression rate is precisely tuned by taking into account this correlation degree and the 
quality of the relay-destination link. In our setup, we have two different receivers. The relay can either 
create a single quantized version of its observation, common to both receivers, or two quantized versions, 
one adapted for each destination. We have chosen the second type of quantization which we call the "bi- 
level compression EF". We note the work by lTT9l where the authors consider a different channel, namely a 
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separated two-way relay channel, and exploit a similar idea, namely using two quantization levels at the relay. 
In the scheme used here, each receiver decodes independently its own message, which is less demanding 
than a joint decoding scheme in terms of information assumptions. As we have already mentioned, the 
relay implements the Wyner-Ziv compression and superposition coding similarly to a broadcast channel. 
The difference with the broadcast channel is that each destination also receives the two direct signals from 
the source nodes. The rate region which can be obtained by using such a coding scheme is given by the 
following theorem proved in Appendix [A] 

Theorem 3.2: For the Gaussian IRC with private messages and bi-level compression EF protocol, any 
rate pair R 2 ) is achievable where 

0j>, (0) subject to the constraints 

AT (1) ^ {\h 11 \ 2 P 1 + \h 21 \ 2 P 2 +N 1 )A-A 2 1 

wz — \h rl \ 2 ^P r 

,.(2) . (\h22\ 2 P2 + \h 12 \ 2 Pi + \h r 2\ 2 viP r +N 2 )A-A 2 2 

wz ~ \h r2 \ 2 v 2 P r ' 

2) dSe ' if ° { \h^P, + \h^ + ^L\^P r + N a ) " ° ( l^PA + I^A+Lp^ + ivJ ' We Mve tlW mteS 
in subject to the constraints 

AT (1) (\hn\ 2 P! + \h2l\ 2 P2 + \h rl \ 2 ^Pr+N^A-Aj 

wz - \h rl \ 2 ^p r 

M (2) . {\h 22 \ 2 P 2 + \h 12 \ 2 Pi+N 2 )A-A 2 2 
wz — \h r2 \ 2 » 2 P r ' 

3) else the rates are given by ((§]), © subject to the constraints 

v (l) (\h 11 \ 2 P 1 + \h 21 \ 2 P 2 + \h rl \ 2 v 2 P r +N 1 )A-A 2 1 

wz — [TvTpI^Tv 

v (2) (\h 22 \ 2 P2 + \h 12 \ 2 P 1 + \h r2 \ 2 iy 1 P r +N 2 )A-A 2 2 

wz — \h r2 \ 2 v 2 P r • 

with Nibl representing the quantization noise corresponding to receiver i, (vi,u 2 ) e [0, l] 2 ,vi + v 2 < 1, the relay 

PA, A = \h lr \ 2 P 1 + \h 2r \ 2 P 2 + N r , Ai = 2TZe{hxih\ r )Pi + 2Ke(h 21 h* 2r )P 2 and A 2 = 2Ke(h 12 h* lr )P 1 + 2Ke(h 22 h* 2r )P 2 . 



January 12, 2011 



DRAFT 



8 



Ri < C 



|/in| 2 F 



V 

/ 



iVl + \h rl \ 2 U 2 Pr + 



\h 2 i\ 2 P2{N r +N, l 



r(l) 



+ 



l^lrP-Pl 



AT j_ ArW _L_ l^2r| 2 P2(|ferl| 2 ^ 2 Pr+A r i) 

iV r -+- 1M WZ -+- T^pT^^J^i^vT 



|/l2r| 2 P 2 +A r r + Afi 1 i 



i?2 C 



No 



|^22| 2 -P2 
Ihial'P^JVr+iV^) 
|fti r | 2 Pi+Ar r +Afi 2 2 



+ 



\h 2r \ 2 P 2 



A r(2) , |fei r | 2 PijV 2 
+ \h 12 \ 2 P l +N 2 



N r 



(6) 



(7) 



/ 



Ri < C 



\ 
( 



Ri C 



No 



\h n \ 2 P 1 



Ni + \h rl \ 2 v 2 Pr 



\h2i\ 2 P2[N r +M, 



r(l) 



\h 2r \ 2 P 2 +N r +N, 



(i) 



|^22| 2 -P2 



\h r2 \ 2 viP r 



2 P 1 (N r +N, 



\h lr \ 2 P!+N r + N, 



(2) 

wz 



N r 



|^lr| 2 -Pl 

ATW I \h2r\ 2 P 2 (\hrl\ 2 V2Pr+N 1 ) 

^ WZ ~f \h2i\ 2 P 2 + \h rl \ 2 v 2 P r +N 1 



+ 



|^2r| 2 -f2 



(8) 



(9) 



The three scenarios emphasized in this theorem correspond to the following situations: l)T>i has the better 
link (in the sense of the theorem) and can decode both the relay message intended for V 2 and its own 
message; 2) this scenario is the dual of scenario 1); 3) in this latter scenario, each destination node sees the 
cooperation signal intended for the other destination node as interference. 

C. Transmission rates for the AF protocol 

In this section, the relay is assumed to implement an analog amplifier which does not introduce any delay 
on the relayed signal. The main features of AF-type protocols are well-known by now (e.g., such relays 
are generally cheap, involve low complexity relay transceivers, and generally induce negligible processing 
delays in contrast with DF and EF-type relaying protocols). The relay merely sends X r = a r Y r where 
a r corresponds to the relay amplification factor/gain. We call the corresponding protocol the zero-delay 
scalar amplify-and-forward (ZDSAF). The type of assumptions we make here fits well to the setting of 
DSL or optical communication networks. In wireless networks, the assumed protocol can be seen as an 
approximation of a scenario with a relay equipped with a power amplifier only. The following theorem 
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provides a region of transmission rates that can be achieved when the transmitters send private messages to 
their respective receivers, the relay implements the ZDSAF protocol and the receivers implement single-user 
decoding. The considered framework is attractive in the sense that an AF-based relay can be added to the 
network without changing the receivers. 

Theorem 3.3 (Transmission rate region for the IRC with ZDSAF): Let R it i e {1,2}, be the transmission 
rate for the source node Si. When ZDSAF is assumed the following region is achievable: 

V? £ {1,2}, R^ F < C[ K h l hri + hii \ 2 P i — — \ (10) 

\Ja r hjrh r i+hji\ pjj^r+a^hri^jfr+l J 

where pi = j^, j = —i, and a r is the relay amplification gain. 

The proof of this result is standard ll20l and will therefore be omitted. Only two points are worth being 
mentioned. First, the proposed region is achieved by using Gaussian codebooks. Second, the choice of the 
value of the amplification gain a r is not always straightforward. In the vast majority of the papers available 
in the literature, a r is chosen in order to saturate the power constraint at the relay (E|X r | 2 = P r ) that 
is: a r = a r = J = J ' \h lr \^p 1 +\h 2r \^p 2 +N r - However, as mentioned in some works EHE3I23E3I. 
this choice can be sub-optimal in the sense of certain performance criteria. The intuitive reason for this is 
that the AF protocol not only amplifies the useful signal but also the received noise. This effect can be 
negligible in certain scenarios for the standard relay channel but can be significant for the IRC. Indeed, 
even if the noise at the relay is negligible, the interference term for user i (i.e., the term hj r Xj, j = —i) 
can be influential. In order to assess the importance of choosing amplification factor a r adequately (i.e., 
to maximize the transmission rate of a given user or the network sum-rate) we derive its best value. The 
proposed derivation differs from ||2T T | ll23l because, here, we consider a different system (an IRC instead 
of a relay channel with no direct link), a specific relaying function (linear relaying functions instead of 
arbitrary functions) and a different performance metric (individual transmission rate and sum-rate instead of 
raw bit error rate ||2T| and mutual information f23J). Our problem is also different from |]24l| since we do 
not consider the optimal clipping threshold in the sense of the end-to-end distortion for frequency division 
relay channels. At last, the main difference with ll22l is that, for the relay channel, the authors discuss the 
choice of the optimal amplification gain in terms of transmission rate for a vector AF protocol having a 
delay of at least one symbol duration; here we focus on a scalar AF protocol with no delay and a different 
system namely the IRC. In this setup, we have found an analytical expression for the best a r in the sense 
of Ri(a r ) for a given user i e {1,2}. We have also noticed that the a r maximizing the network sum-rate 
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has to be computed numerically in general. The corresponding analytical result is stated in the following 
theorem. 

Theorem 3.4: [Optimal amplification gain for the ZDSAF in the IRC] The transmission rate of user i, 
Ri(a r ), as a function of a T G [0,a r ] can have several critical points which are the real solutions, denoted 
by cfl and cf], to the following second degree equation: 

a 2 [\mi\ 2 Re(piq*) - (\pi\ 2 + s^Re^n*)] + a r [|m;| 2 (|g;| 2 + 1) - K| 2 (|p;| 2 + S;)] 
+ (|%| 2 + l)Re(m^*) - \ni\ 2 Re{ Pi q*) = 

where rrii = h ir h riy /p~, rii = h iiy /p~, p { = h jr h riy /p], q { = h jiy /pj, = \h ri \ 2 , i G {1,2} and j = -i. 
Thus, depending on the channel parameters, the optimal amplification gain a* = arg max Ri(a r ) takes 

a r G[0,a r ] 

one value in the set a* G {0, a n c^j, c^j}. If, additionally, the channel gains are reals then the two critical 
points write as: c$ = and cf] = - . 

The proof of this result is provided in Appendix [Bj Of course, in practice, if the receive signal-to-noise 
plus interference ratio (viewed from a given user) at the relay is low, choosing the amplification factor a r 
adequately does not solve the problem. It is well known that in real systems, a more efficient way to combat 
noise is to implement error correcting codes. This is one of the reasons why DF is also an important relaying 
protocol, especially for digital relay transceivers for which AF cannot be implemented in its standard form 
(see e.g., |[24l for more details). 



D. Time-Sharing 

In terms of achievable Shannon rates, distributed channels differ from their centralized counterpart. Indeed, 
rate regions are not necessarily convex since the time-sharing argument can be invalid (if no synchronization 
is possible). Similarly, depending on the channel gains, the achievable rate for a given transmitter can be non- 
concave with respect to its power allocation policy. This happens if the transmitters cannot be coordinated 
(distributed channels). 

Assuming that the users can be coordinated, we discuss here a standard time-sharing procedure similarly 
to the approach in E51 for the frequency-division relay channel. More specifically, we assume that user 1 
decides to transmit only during a fraction ol\ of the time using the power ^ and user 2 transmits only with 
a fraction a 2 percent of the time using the power ||. 

The achievable rate-region with coordinated time-sharing, irrespective of the relay protocol, is: 
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Vz G {1,2}, Rj s < o^i* ( — , ) + Oi^-iJi (-, ^ ) , (12) 

where j = —i, (a i: a,) 2 G [0, l] 2 , (/?,, /i,-) 2 G [0, l] 2 such that /^a^ = fiiOt\. The rate Oj represents 

the achievable rate of user i (depends on the relay protocol and was provided in the previous subsections) 
when user j doesn't transmit and user i transmits with power — , E4 ( — , — ) is the achievable rate when 

P P 

user % transmits with power and user j transmits with power In order to achieve this rate region, 
the users have to be coordinated. This means that they have to know at each instant if the other user is 
transmitting or not. User % also has to know the parameters on and cij. The parameter f3j G [0,1] represents 
the fraction of time when user j interferes with user i. Considering the time when both users transmit with 
non-zero power, we obtain the condition: /3ia 2 = Ajai. 

IV. Power allocation games in multi-band IRCs and Nash equilibrium analysis 

In the previous section, we have considered the system model presented in Sec. [II] for Q = 1. Here, 
we consider multi-band IRCs for which Q > 2. As communications interfere on each band, choosing the 
power allocation policy at a given transmitter is not a simple optimization problem. Indeed, this choice 
depends on what the other transmitter does. Each transmitter is assumed to optimize its transmission rate 
in a selfish manner by allocating its transmit power Pj between Q sub-channels and knowing that the other 
transmitters want to do the same. This interaction can be modeled as a strategic form non-cooperative 
game, Q = (/C, (Aiji&c, (i*i)i e /c), where: (i) the players of the game are the two information sources or 
transmitters and /C = {1,2} is used to refer to the set of players; (ii) the strategy of transmitter % consists 

9^ < 1 > where t - represents 



in choosing ^ = (B\ \ . . . , 6>j g) ) in its strategy set At = < ^ G [0, 1] Q 

the fraction of power used in band (q); (iii) itj(-) is the utility function of user % G {1,2} or its achievable 
rate depending on the relaying protocol. From now on, we will call state of the network the (concatenated) 
vector of power fractions that the transmitters allocate to the IRCs i.e., 9 = {9^,9^). An important issue 
is to determine whether there exist some outcomes to this conflicting situation. A natural solution concept 
in non-cooperative games is the Nash equilibrium [26|. In distributed networks, the existence of a stable 
operating state of the system is a desirable feature. In this respect, the NE is a stable state from which the 
users do not have any incentive to unilaterally deviate (otherwise they would lose in terms of utility). The 
mathematical definition is the following. 
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Definition 4.1 : [Nash equilibrium] The state J is a pure NE of the strategic form game Q if 

Mi e JC,W, e A, uMliU) > 

In this section, we mainly focus on the problem of existence of such a solution, which is the first step 
towards equilibria characterization in IRCs. The problems of equilibrium uniqueness, selection, convergence, 
and efficiency are therefore left as natural extensions of the work reported here. 

A. Equilibrium existence analysis for the DF protocol 

As explained in Sec. IIII-AI the signals transmitted by Si and <5>2 m band (q) have the following form: 
X[ q) = X$ + a T \ q) 6i \q) X^ where the signals X$ and are Gaussian and independent. At the relay 

y v i p r ' 

71, the transmitted signal writes as: X^ = X*fJ + X^. For a given allocation policy 6 { = (o^, Q^J, 
the source-destination pair (Si,Vi) achieves the transmission rate Y^a=i DF wnere 



i^ DF = min{i?S 9 J' DF ,i?g' DF } 
R^ F = min{i?g' DF ,i?g' DF } 



(13) 



and Rfi° F , R$ BF ,R^ DF ,R%% DF are given in (fl4l) and (v( q \ r[ 9 \ r 2 ) is a given triple of parameters in 
[0,l]\r[ q) +4 q) <l. 
The achievable transmission rate of user i is given by: 

u^U^R^ie^^j). (15) 

9=1 

We suppose that the game is played once (one-shot or static game), the users are rational (each selfish player 
does what is best for itself), rationality is common knowledge, and the game is with complete information 
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that is, every player knows the triplet Q DF = (JC, (A)ie£! (uf F ) ie )c)- Although this setup might seem to be 
demanding in terms of CSI at the source nodes, it turns out that the equilibria predicted in such a framework 
can be effectively observed in more realistic frameworks where one player observes the strategy played by 
the other player and reacts accordingly by maximizing his utility, the other player observes this and updates 
its strategy and so on. We will come back to this later on. The existence theorem for the DF protocol is 
given hereunder. 

Theorem 4.2: [Existence of an NE for the DF protocol] 
If the channel gains satisfy the condition 1Ze(h\f h^*) > 0, for all i G {1, 2} and q G {1, . 

defined by G DF = (/C, (AW, (u? F (0i,£_i)W) with K = {1,2} and A = l^E [0, 1] ( 



Q "| 

£«?><! , 

9=1 J 



has always at least one pure NE. 

Proof: The proof is based on Theorem 1 of Il2~71l . The latter theorem states that in a game with a finite 
number of players, if for every player 1) the strategy set is convex and compact, 2) its utility is continuous 
in the vector of strategies and 3) concave in its own strategy, then the existence of at least one pure NE 
is guaranteed. In our setup checking that conditions 1) and 2) are met is straightforward. The condition 
Tle(hif h^i*) > is a sufficient condition that ensures the concavity of Rf F w.r.t. 6^. The intuition behind 
this condition is that the superposition of the two signals carrying useful information for user % (i.e., signal 
from Si and 1Z) has to be constructive w.r.t. the amplitude of the resulting signal. As the sum of concave 
functions is a concave function, the min of two concave functions is a concave function (see ll28l for more 



details on operations preserving concavity), and B!fj is a concave function of 6 { , it follows that 3) is also 
met, which concludes the proof. ■ 
Theorem indicates, in particular, that for the pathloss model where the channel gains are positive real 
scalars (i.e., hij > 0, G {l,2,r} 2 ) there always exists an equilibrium. As a consequence, if some 

relays are added in the network, the transmitters will adapt their PA policies accordingly and, whatever 
the locations of the relays, an equilibrium will be observed. This is a nice property for the system under 
investigation. As the PA game with DF is concave it is tempting to try to verify whether the sufficient 
condition for uniqueness of ll27l is met here. It turns out that the diagonally strict concavity condition of 
ll27l is not trivial to be checked. Additionally, it is possible that the game has several equilibria as it is 
proven to be the case for the AF protocol. 
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5. Equilibrium existence analysis for the EF protocol 

In this section, we make the same assumptions as in Sec. IIV-AI concerning the reception schemes and 
PA policies at the relays: we assume that each node 1Z, T>i and V 2 implements single-user decoding and 
the PA policy at each relay i.e., y_ = (V 1 * 1 , u^) is fixed. Each relay now implements the EF protocol. 
Under this assumption, the utility for user % e {1,2} can be expressed as follows 



2? (?),EF 
q=l 

where are given in (O G [0, 1], A® = \h { $\^ q) P^h®]^ P 2 +N^ , A[ q) = h^h^'*8[ q) Pi+ 

h { 2 q }h { 2 f*8 { 2 q) P 2 and A { 2 q) = h$ h$'* 8 {q) P x + h 2 q J h 2 f* 8 (q) P 2 . What is interesting with this EF protocol is 
that, here again, one can prove that the utility is concave for every user. This is the purpose of the following 
theorem. 

Theorem 4.3: [Existence of an NE for the bi-level compression EF protocol] The game defined by Q EF = 
(JC,(A i ) ie K,(y'f F (li,0- i ))ie!c) with K = {1,2} and A = E [0, 1]« 
one pure NE. 

The proof is similar to the proof of Theorem IIV-AI To be able to apply Theorem 1 of Rosen [27], we 
have to prove that the utility uf F is concave w.r.t. 8^. The problem is less simple than for DF because the 
compression noise n!~ i which appears in the denominator of the capacity function in Eq. (fTTT) depends on 
the strategy 8 { of transmitter i. It turns out that it is still possible to prove the desired result as shown in 
Appendix O 
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C. Equilibrium analysis for the AF protocol 

Here, we assume that the relay implements the ZDSAF protocol, which has already been described in Sec. 
IIII-Cl One of the nice features of the (analog) ZDSAF protocol is that relays are very easy to be deployed 
since they can be used without any change on the existing (non-cooperative) communication system. The 
amplification factor/gain for the relay on band (g) will be denoted by of 1 . Here, we consider the most 
common choice for the amplification factor that it, the one that exploits all the transmit power available on 
each band. The achievable transmission rate is given by 



q=l 



(19) 



where R^' AF is the rate user i obtains by using band (q) when the ZDSAF protocol is used by the relay 
TZ. After Sec. IIII-Cl the latter quantity is: 



\/ie{l,2},R, 



where = $\e[ 9 \ 9?) 



(<?),AF 



c 



\ 



|4«> fcirftrl+hji | 3 pj9 («)^! 7+ ( («)) a | fc (5)| a 
\ i i 



(20) 



Pr 



Pl + |/l2r| 2 P2 + iV r 



and pf' = -^y. Without loss of generality and for the 



sake of clarity we will assume in Sec. HV-Cl that V(i, q) G {1, 2, r} x {1, . . . , Q}, N> 9) — N, P r w = P r and 
we introduce the quantities = In this setup the following existence theorem can be proven. 



Theorem 4.4: [Existence of an NE for ZDSAF] If any of the following conditions are met: i) 
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1 %r ri 
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and 



(negligible direct links), ii) 
(negligible relay links), Hi) aft = Ar G [0, Or (1,1)] (constant amplifi- 
cation gain), there exists at least one pure NE in the PA game Q AF . 

The proof is similar to the proof of Theorem IIV-AI The sufficient conditions ensure the concavity of the 



function R^' AF w.r.t. Of ■ For the first case i) where the direct links between the sources and destinations 
are negligible (e.g., in the wired DSL setting these links are missing and the transmission is only possible 
using the relay nodes), the achievable rates are given Vz G {1, 2}, in (T2TI) and it can be proven that Rr"' AF is 



j(9) 
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concave w.r.t. #| . The other two cases are easier to prove since the amplification gain is either constant or 
not taken into account and the rate R^' AF is a composition of a logarithmic function and a linear function 
of 9^ and thus concave. 

The determination of NE and the convergence issue to one of the NE are far from being trivial in this case. 
For example, potential games ||29l and supermodular games 113011 are known to have attractive convergence 
properties. It can be checked that, the PA game under investigation is neither a potential nor a supermodular 
game in general. To be more precise, it is a potential game for a set of channel gains which corresponds 
to a scenario with probability zero (e.g., the parallel multiple access channel). The authors of OTl studied 
supermodular games for the interference channel with K = 2, Q = 3, assuming that only one band is shared 
by the users (IC) while the other bands are private (one interference-free band for each user). Therefore, each 
user allocates its power between two bands. Their strategies are designed such that the game has strategic 
complementarities. However, as stated in OTl . this design trick does not work for more than two players or 
if the users can access more than two frequency bands. In conclusion, general convergence results seem to 
require more advanced tools and further investigations. 

Special case study 

As we have just mentioned, the uniqueness/convergence/efficiency analysis of NE for the DF and EF 
protocols requires a separate work to be treated properly. However, it is possible to obtain relatively easy 
some interesting results in a special case of the AF protocol. The reason for analyzing this special case 
is threefold: a) it corresponds to a possible scenario in wired communication networks; b) it allows us to 
introduce some game-theoretic concepts that can be used to treat more general cases and possibly the DF and 
EF protocols; c) it allows us to have more insights on the problem with a more general choice for af \ The 
special case under investigation is as follows: Q = 2 and Vg £ {1, 2}, ai = A^ £ [0, a r (l, 1)] are constant 
w.r.t. 9. We observe that the strategy set of user i are scalar spaces 8i £ [0, 1] because we can consider 
0\ = 9i and of' = 8^. For the sake of clarity, we denote by = h$ and = Note that the case 
af^ = A^ can also be seen as an interference channel for which there is an additional degree of freedom 
on each band. The choice Q = 2 is totally relevant in scenarios where the spectrum is divided in two bands, 
one shared band where communications interfere and one protected band where they do not (see e.g., 113210 . 
The choice a?^ = const, has the advantage of being mathematically simple and allows us to initialize the 
uniqueness/convergence analysis. Moreover, it corresponds to a suitable model for an analog repeater in 
the linear regime in wired networks or, more generally, to a power amplifier for which neither automatic 
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gain control is available nor received power estimation mechanism. By making these two assumptions, it 
is possible to determine exactly the number of Nash equilibria through the notion of best response (BR) 
functions. The BR of player i to player j is defined by BR^fy) = argmax«j(0j, 9j). In general, it is a 
correspondence but in our case it is just a function. The equilibrium points are the intersection points of 
the BRs of the two players. In this case, using the Lagrangian functions to impose the power constraint, it 
can be checked that: 

' Fjfij) if < FiiOj) < 1 

1 if Fitfj) > 1 (22) 
otherwise 



BRj(0j, 



where j = -i, F^Oj) = is an affine function of 6j for (ij) £ {(1,2), (2,1)}, c tt = 2\A ( r 1) h ri h ir + 

ha\ \A r QriQir ~\~ 9ii\ Pit = \Ar h r j}lj, r -\- \A r 9ri9jr~\~9ji\ Pj ~\~\A r fl r ihj r + hji\ \A r QriQir ~\~ 9ii\ Pj* 

di = \Ar ^ h r ihi r -\-hii\ 2 [\ Ar ^ g r i9iT~\~ 9ii\ 2 Pi~^~\^r ^ 9ri9jr~\~ 9ji\ 2 Pj + ^ \9ri \ 2 + 1] — | Ar ' ^ 9ri9ir~\~ 9ii\' ''(A r ^ \h r i\ 2 + 

1). By studying the intersection points between BRi and BR 2 , one can prove the following theorem (the 
proof is provided in Appendix iDl). 

Theorem 4.5 (Number of Nash equilibria for ZDSAF): For the game Q AF with fixed amplification gains 
at the relays, (i.e., = 0), there can be a unique NE, two NE, three NE or an infinite number of NE, 
depending on the channel parameters (i.e., hij, gij, p^ Ai q \ G {l,2,r} 2 , q £ {1,2}. 

Notice that, if A r = 0, we obtain the complete characterization of the NE set for the two-users two- 
channels parallel interference channel. In the proof in Appendix O we give explicit expressions of the 
possible NE in function of the system parameters (i.e., the amplification gain A r and the channel gains). If 
the channel gains are the realizations of continuous random variables, it is easy to prove that the probability 
of observing the necessary conditions on the channel gains for having two NEs or an infinite number of 
NEs is zero. Said otherwise, considering the pathloss model and arbitrary nodes positioning, there will be, 
with probability one, either one or three NE, depending on the channel gains. When the channel gains are 
such that the NE is unique, the unique NE can be shown to be: 

flNE = Q* = ( C 22Q?1 ~ C 12 d 2 Ciirf 2 ~ C 2 idi \ 



.CllC 22 - C12C21 C11C22 - C12C21, 

When there are three NE, it seems a priori impossible to predict the NE that will be effectively observed in 
the one-shot game. In fact, in practice, in a context of adaptive/cognitive transmitters (note that what can be 
adapted is also the PA policy chosen by the designer/owner of the transmitter), it is possible to predict the 



January 12, 2011 



DRAFT 



18 



equilibrium of the network. First, in general, there is no reason why the sources should start transmitting 
at the same time. Thus, one transmitter, say i, will be alone and using a certain PA policy. The transmitter 
coming after, namely <S_j, will sense/measure/probe its environment and play its BR to what it observes. As 
a consequence, user i will move to a new policy, maximizing its utility to what transmitter — i has played 
and so on. The key question is: does this procedure converge? This procedure is guaranteed to converge to 
one of the NE and a detailed discussion about the asymptotic stability of the NE can be found in Appendix 
151 The arguments for proving this have been used for the first time in ll33l where the "Cournot duopoly" was 
introduced. In 11331 . the BRs of each player is purely affine, which leads in this case to a unique equilibrium. 
The corresponding iterative procedure is called the Cournot tatonnement process in 11341 . In the case with 
three NE, the effectively observed NE can be predicted by knowing the initial network state that is, the 
PA policy played by the first transmitting player (see Sec. |V]). To implement such an iterative procedure, it 
can be checked [1] that the transmitters need to know less network parameters than in the original game 
where the amplification factor saturates the constraint. In fact, the needed parameters can be acquired by 
realistic sensing/probing techniques or feedback mechanisms based on standard estimation procedures. As 
a comment, note that in the (modern) literature of decentralized or distributed communications networks 
where the optimal PA policy of a transmitter is to water-fill, the mentioned iterative procedure is called 
iterative water-filling. 

D. Equilibrium analysis for the Time-Sharing scheme 

In the previous subsections, we have given sufficient conditions that ensure the existence of the Nash 
equilibrium. Our approach is based on the concave games studied in E71 and consists in finding the sufficient 
conditions that ensure the concavity of the transmission achievable rates. We have seen that, assuming 
ZDSAF or DF relaying protocols, the achievable rates are not necessarily concave. 

Assuming that the transmitters can be coordinated, and by using the time-sharing scheme similarly to 
Subsec. IIII-Dl the achievable transmission rate of user i is given by: 

^fiJ = E^' TS (fiS fl) ^i ) ) (24) 

q=l 

where i?^' TS is the rate user i obtains by using band (q) and time-sharing technique. After Sec. IIII-Dl the 
latter quantity is: 

Vi € {1, 2}, fl<* TS = a' V*,!" o) + af'pf'^ , (25) 
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where j = -i, (a?\af) e [0, l] 2 , ($ q) , fif) e [0, l] 2 such that p[ q) a$ = ^ 2 q) af\ These parameters are 

( ) ( ffe^ p i ( \ ( p o ^ p m \ 

fixed and chosen such that the achievable rates are maximized. The rates Rf' ( -» — i , I , Ry 1 ( -» — i , — - I 

1 I O, ' I' « I OLi OLj J 

represent the achievable rates in band (q) when time-sharing is used. These rates depend on the relaying 
protocol and are given by Eq. (fl~3l) for DF and by Eq. (1201) for ZDSAF. Notice that, when EF is assumed, 
the rates are always concave irrespective of the channel gains and time-sharing techniques do not change 
the achievable rate-region. 

Theorem 4.6: [Existence of an NE for TS] There always exists at least one pure NE in the PA game Q TS , 
regardless of the used relaying scheme and the values of the channel gains. 

If the users are coordinated (i.e., each user is aware of the moments where the other user is transmitting 
or not) then their achievable rates PSf 1 ^ are always concave w.r.t. 6\ q . This implies directly that 11271 . 
irrespective of the relaying technique and of the channel gains, the existence of an NE will be guaranteed. 

In the particular case where either F r (g) = or h {q ) = 0, for all q G {1, . . . , Q} and i G {1, 2}, the parallel 
IRC reduces to the parallel interference channel [6]. The time-sharing scheme is useless since the achievable 
rates are already concave and a, = 1, $ = 1 are optimal. Therefore, Theorem 14.61 guarantees the existence 
of the NE in this case and is consistent with the known results in J6J. 

V. Simulation results 

Single-band IRCs: AF vs DF vs EF. Here, we assume Q — 1 and a path loss exponent of 2 that is, 
\hij\ = 0^) 2 f° r (hj) ^ {l)2,r} 2 where d Q = 5 m is a reference distance and 7 = 2 is the path loss 
exponent. The nodes <Si, S 2 , T>\, V 2 are assumed to be located in a plane. The positions of the nodes will be 
indicated on each figure and are characterized by the distance between them which are chosen as follows: 
d' u = 11.5 m, d' 22 = 10 m, d' 12 = 11 m and d' 21 = 14 m. As for the relay, to avoid any divergence for the 
path loss in dy = 0, we assume that it is located a hight e = 0.1 m from this plane i.e., the relay location 



is given by the (x r , y r , z r ) where z r is fixed and equals 0.1m; thus dij = yd'£j + e 2 for 2 = r or j = r and 
% 7^ j. The noise levels at the receiver nodes are assumed to be normalized (Ni = N 2 = N r = 1). In terms 
of transmit power we analyze two cases: a symmetric case where Pi = P 2 = 10 (normalized power) and an 
asymmetric one where P\ = 3 and P 2 = 10. The relay transmit power is fixed: P r = 10. For the symmetric 
scenario, Fig. |2] represents the regions of the plane f^^H G [ — 4, +4] x [—3, +4] (corresponding to the 
possible relay positions) where a certain protocol performs better than the two others in terms of system 
sum-rate. These regions are in agreement with what is generally observed for the standard relay channel. 
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This type of information is useful, for example, when the relay has to be located in specific places because 
of different practical constraints and one has to choose the best protocol. Fig. fallows one to better quantify 
the differences in terms of sum-rate between the AF, DF and bi-level EF protocols since it represents the 
sum-rate versus x r for a given y r = 0.5d . The discontinuity observed stems from the fact that for the 
bi-level EF protocol there is a frontier delineating the scenarios where one receiver is better than the other 
and can therefore suppress the interference of the relay (as explained in Sec. IIII-B|) . 

Number of Nash equilibria for the AF protocol. First, we show that in the PA game with ZDSAF, one can 
have three possible Nash equilibria. For a given typical scenario composed of an IC in parallel with an IRC 
(Q = 2) and pi = 1, p 2 = 3, p r = 2 and the channel gains (gn, g\ 2 , g 2 i, #22) = (2.76,5.64,-3.55,-1.61), 
(hii,h 12 ,h 21 ,h 22 ) = (14.15,3.4,0,1.38) and (h lr , h 2r , h rl , h r2 ) = (-3.1, 2.22, -3.12, 1.16), we plot 
the best response functions in Fig. HI We see that there are three intersection points and therefore three Nash 
equilibria. As explained in Sec. IIV-CI the effectively observed NE in a one-shot game is not predictable 
without any additional assumptions. However, the Cournot tatonnement procedure converges towards a given 
NE which can be predicted from the sole knowledge of the starting point of the game, namely 9\ or 9®. 

Stackelberg formulation. We have mentioned that a strong motivation for studying IRCs is to be able to 
introduce relays in a network with non-coordinated and interfering pairs of terminals. For example, relays 
could be introduced by an operator aiming at improving the performance of the communications of his 
customers. In such a scenario, the operator acts as a player and more precisely as a game leader in the 
sense of ll35l . In ll35l . the author introduced what is called nowadays a Stackelberg game. This type of 
hierarchical games comprises one leader which plays in the first step of the game and several players (the 
followers) which observe the leader's strategy and choose their actions accordingly. In our context, the 
game leader is the operator/engineer /relay who chooses the parameters of the relays. The followers are the 
adaptive/cognitive transmitters that adapt their PA policy to what they observe. In the preceding sections we 
have mentioned some of these parameters: the location of each relay; in the case of AF, the amplification gain 
of each relay; in the case of DF and EF, the power allocation policy between the two cooperative signals at 
each relay i.e., the parameter v^ q \ Therefore, the relay can be thought of as a player who maximizes its own 
utility. This utility can be either the individual utility of a given transmitter (picture one WiFi subscriber 
wanting to increase his downlink throughput by locating his cellular phone somewhere in his apartment 
while his neighbor can also exploit the same spectral resources) or the network sum-rate (in the case of an 
operator). In the latter case, the operator possesses some degrees of freedom to improve the efficiency of 
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the equilibrium. In the remaining part of this section, we focus on the Stackelberg formulation where the 
strategy of the leader is respectively the relay amplification factor, position and power allocation between the 
cooperative signals. The system considered is composed of an IRC in parallel with an interference channel 
(IC) [36]. All the simulations provided are obtained by applying the Cournot tatonnement procedure. The 
simulation setup is as follows. The source and destination nodes are located in fixed locations in the region 
[— L, L] 2 of a plane, with L = 10 m, such that the relative distances between the nodes are: dn = 6.52m, 
di2 = 8.32m, di\ = 6.64m, d 22 = 6.73m. We assume a path loss model for the channel gains \g%j\, \hu\. 

7 (1) 7 (2) 

For the path loss model we take \hij\ = O^j 2 an d \g%j\ = {^J 2 f° r (hj) e {l>2,r} 2 where 
do = 1 m is a reference distance. The relay is at e = 0.5 m from the plane. We will also assume that 

N (i) = N (2) = N ^ . G |^ 2 ^ ^(l) wm be denoted by Nr and also p (i) = P ^ A m = A ^ a a) = a ^ 

= a r , = v. 

Optimal relay amplification gain for the AF protocol. First we consider the ZDSAF relaying scheme 
assuming a fixed amplification gain a r = A r (Sec. IIV-CI ). We want to analyze the influence of the value of 
the amplification factor, A r £ [0, a r (l, 1)], on the achievable network sum-rate at the NE. This is what Fig. 
[5] shows for the following scenario: e = 0.5 m, Pi = 20 dBm, P 2 = 23 dBm, P r = 22 dBm, N\ = 10 dBm, 
N 2 = 9 dBm, N r = 7 dBm, 7 (1) = 7 (2) = 2. We observe that the optimal value is A* = 0.05 and is not 
equal to the one saturating the relay power constraint a r (l, 1) = 0.17. This result illustrates for the sum-rate 
what we have proved analytically for the individual rate of a given user (see Sec. IIII-CI) . Note that the gap 
between the optimal choice for a r and the choice saturating the power constraint is not that large and in 
fact other simulation results have shown is is generally of this order and even smaller typically. 

Optimal relay location for the AF protocol. Now, we consider the ZDSAF when the full power regime 
is assumed at the relay, o r — (i r {B\ , 62) (Sec. IIV-CI) and study the relay location problem. Fig. [6] represents 
the achievable network sum-rate as a function of the relay position (xiz,yn) G [—L,L] 2 for the scenario: 
Pi = 20 dBm, P 2 = 17 dBm, P r = 22 dBm, iVi = 10 dBm, N 2 = 9 dBm, N r = 7 dBm, 7W =2.5 
and 7^) = 2. We observe that there are two local maximum that actually correspond to the points that 
maximize the individual achievable rates. Many simulation results have confirmed that, when the source 
nodes are sufficiently far away from each other, maximizing the individual rate of either user at the NE 
amounts to locating the relay on one of the the segment between <S, and T>i. This interesting and quite 
generic observation can be explained as follows. For this purpose, consider Fig. [6] which is a temperature 
image representing the values of 0\ and 9 2 for different relay positions in [— L, L} 2 . The region where 
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(#i,6> 2 ) = (1,0) (resp. (9i,6 2 ) = (0,1)) is the region around Si (resp. <S 2 ). We see that the intersection 
between these regions corresponds to a small area. This quite general observation shows that the selfish 
behavior of the transmitters leads to self-regulating the interference in the network. Said otherwise, a selfish 
transmitter will not use at all a far away relay but leaves it to the other transmitter. Thus, when one transmitter 
uses the relay, it is often alone and sees no interference. In these conditions, by considering the path loss 
effects it can be proved that the optimal relay position is on the segment between the considered source 
and destination nodes. This also explains why the position that maximizes the network sum-rate lies also 
on one of the segments from Si to 

Optimal relay power allocation at the relay for DF and EE For the DF protocol, we fix the cooperation 
degrees T\ — and r 2 = 0. In Fig. [7J we plot the achievable sum-rate at the equilibrium as a function of 
the relay power allocation policy is v £ [0, 1] (with the convention v = v^) for the scenario: x-ji = m, 
y n = m, Pi = 22 dBm, P 2 = 17 dBm, P r = 23 dBm, Ni = 7 dBm, N 2 = 9 dBm, N r = dBm, 
= 2.5 and 7^ = 2. We observe that, for both protocols, the optimal power allocation v* = 1, meaning 
that the relay allocates all its available power to the better receiver, T>i. In this case, the relay is in very good 
conditions and can therefore reliably decode the source messages. This explains why DF outperforms EF 
which is in agreement with the observations we have made in Sec. Hill We have observed that, in general, 
the network sum-rate is not concave w.r.t. v £ [0, 1] and that the optimal power allocation lies on the borders 
v* £ {0, 1} for both relaying protocols. In Fig. |TJ we also see that the fair PA policy that is, v — \ can lead 
to a relatively significant performance loss. 

VI. Conclusion 

The complete study of PA games in IRCs is a wide problem and we do not claim to fully characterize it 
here. One of the main objectives in this paper has been to know whether there exist some stable outcomes 
to the conflicting situation where two transmitters selfishly allocate their power between different sub- 
channels in multi-band interference relay channels in order to maximize their individual transmission rate. 
Our approach has been to consider transmission rates achievable in a decentralized framework where relays 
can be deployed with minor or even with no changes for the already existing receivers. For the three types of 
protocols considered, we have proved that the utility of the transmitters is a concave function of the individual 
strategy, which ensures the existence of Nash equilibria in the power allocation game after Rosen ll27l . In a 
special case of the AF protocol, we have fully characterized the number of NE and the convergence problem 
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of Cournot-type or iterative water-filling procedures to an NE. Although we have limited the scope of the 
paper, we have seen that studying IRCs deeply requires further investigations. Many interesting questions 
which can be considered as natural extensions of this work have arisen. Considering more efficient coding- 
decoding schemes and relaying protocols such as those of |[T5l and related works, is it possible to prove 
that the utilities are still concave functions? For these schemes and those considered in this paper, it is also 
important to fully determine the number of Nash equilibria and derive convergent iterative distributed power 
allocation algorithms. We have also seen that several power allocation games come into play and need to be 
studied when considering DF, EF and AF-type protocols: for allocating transmit power between the different 
bands at the sources, for choosing the cooperation degree at the sources, for allocating the power between 
the cooperation signals at the relay, for allocating the transmit power over time. Furthermore, a new agent 
can come into play (the relay) and several Stackelberg formulations can be used to improve the efficiency 
of the equilibria. 
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Appendix A 

Proof of Theorem 13 . 21 (achievable transmission rates for IRCs with the EF protocol) 

In order to prove that the transmission rate region of Theorem 13.21 is achievable for Gaussian IRCs, 
we use a quite common approach ll2~0ll for proving coding theorems: we first prove that it is achievable 
for discrete input discrete output channels and obtain the Gaussian case from standard quantization and 
continuity arguments ll20l . and a proper choice of coding auxiliary variables. 

Definitions and notations 
We denote by Ae (X) the weakly e-typical set for the random variable X. If X is a discrete variable, 
X e X, then || X\\ denotes the cardinality of the finite set X. We use x n to indicate the vector (xi, ■ ■ ■ , x n ). 

Definition A. 1 : A two-user discrete memoryless interference relay channel (DMIRC) without feedback 
consists of three input alphabets X\, X 2 and X r , and three output alphabets 3^i, 3^2 and y r , and a probability 

n 

transition function that satisfies p(yi,V2,y? I x ii x 2i x r) = \\p{Vi,kiy2,k-,yr,k I x \,ki x 2,k, x r,k) fo r some 
n G N*. 

Definition A. 2: A (2 nRl , 2 nFt2 , nj -code for the DMIRC with private messages consists of two sets of 
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integers Wi = {l, 2 nRl } and W 2 = {l, 2 nR2 }, two encoders: /< : -> a se? of relay 

functions {f r ,k}k=i sucn tnat x r,k = fr,k (y r ,i,y r ,2, ■■■,y r ,k-i) , 1 < k < n and two decoding functions 
9% '■ y°i —> i G {1, 2}. 77ze source node Si intends to transmit W it the private message, to the receiver 
node T>i. 

Definition A. 3: The average probability of error is defined as the probability that the decoded message 
pair differs from the transmitted message pair; that is, P^ = Pr [gx (F") ^ W\ or g 2 (Y 2 n ) 7^ W2 \ {W\, W 2 ) 
where (Wi,W 2 ) is assumed to be uniformly distributed over VVi x W 2 . We also define the the average 
probability of error for each receiver as P^ = Pr [gi (YD 7^ W\ I Wj\ . We have < max jpjf\ P^ j < 
pj n) < pj£> + P e ( "l //ence P e (n) -»■ implies that both P^ -> and P e 2 -> 0, an J conversely. 

Definition A.4: A rate pair (Ri,R 2 ) is said to be achievable for the IRC if there exists a sequence of 
(2 nR \2 nR2 ,n) codes with P e (n) ->■ as n -> 00. 

Overview of coding strategy 
At the end of the block k, the relay constructs two estimations y^k) and y^ik) of its observation y^(i) 
that intends to transmit to the receivers T>\ and V 2 to help them resolve the uncertainty on and w 2 .k 
respectively at the end of the block k + 1. 

Details of the coding strategy 
Codebook generation 

i Generate 2 nRl i.i.d. codewords x?(wi) ~ Ylk=i P where i«j G {l, . . . ,2 nRi j, i G {1,2}. 

ii Generate 2 nR ° ' i.i.d. codewords w" ~ Ilfc=iP ( u i,k)- Label these u™(si), s\ G Jl, . . . , 2 n,R ° ' 

iii Generate 2 rtR o i.i.d. codewords u 2 ~ Ilfc=i P ( u 2,k)- Label these m 2 (s 2 ), $i G si,..., 2 nR ° 

iv For each pair (it"(si), m 2 (S2)), choose a sequence x" where x" ~ p (x"|u"(si), u 2 (s 2 )) = 

llfe=lP ( x r,fe|«l,fe(si), M 2 ,fc(s 2 )). 

v For each u"(si), generate 2 nRl conditionally i.i.d. codewords y^ ~ Ilfc=i P (yrik\ui,k(si)) and 
label them ^(zjjsi), z 1 G |l, . . . , 2 n ^ 1 |. For each pair (ux,y rl ) G Ui x the conditional 

probability jp(# r i|wi) is defined as p{y r i\ui) = ^ p(^i)p(a;2)p(yi,y2,?/rki,a;2,a; r )p(y r i|y I .,ui). 

x 1 ,x 2 ,yi,y-2,yr 

vi For each ii 2 (s 2 ), generate 2 nR2 conditionally i.i.d. codewords y™ 2 ~ rifc=iP (yr2k\u2,k(s 2 )) and label 
them yr 2 (z 2 \s 2 ), z 2 G |l,2 n ^ 2 |. For each triplet (u 2l y r i) G U 2 x 5ri 5 the conditional probability 
K#r2|M 2 ) is defined as p(y r2 \u 2 ) = p (xi)p (x 2 )p {yx, y 2 , y r \xi, x 2 ) V (jtall/rj ^2)- 



xi,x2,yi,y2,yr 



vii Randomly partition the message set jl, 2, . . . , 2 n ^| into 2 n " R o 1) sets js^, S^, . . . , ^^(i) | b Y 
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(2) 



independently and uniformly assigning each message in 1 1 , . . . , 2 nRl j to an index in 1 1 , . . . , 2 nR o ) j . 

viii Also, randomly partition the message set |l, 2, . . . , 2 n ^ 2 | into 2 nR o 2> sets jsf^, S {2) r(2) j 

by independently and uniformly assigning each message in 1 1 , . . . , 2 n ^ 2 j to an index in 1 1, ... , 2 n - R o 
Encoding procedure Let w ljk and w 2 ,k be the messages to be send in block k. Si and S 2 respectively 
transmit the codewords x™(wi jk ) and a^C^fc)- We assume that ^(zi^-ilsi^.i), y^(k — 1)) G 

Ae"'' and zi^-i G and also that 

(^(s2,fc-i), ^(^fc-iN.fc-i), £/" — 1)) G with Z2,fe-i £ S^ k . Then the relay transmits the codeword 

<Ol,fc,S 2 ,fc)- 

Decoding procedure In what follows, we will only detail the decoding procedure at the receiver node T>\ 
(at T>2 the decoding is analogous). At the end of block k: 

i The receiver node T>\ estimates s l fc = s\ if and only if there exists a unique sequence u" (si) that 
is jointly typical with Hi(k). We have S\ = s x fc with arbitrarily low probability of error if n is 



sufficiently large and Rq < I{U\\ Yj 



ii Next, the receiver node Vi constructs a set L\ (y™ (k — 1)) of indexes Z\ such that 

(«i (si,fe-i) iVri (^i|si,fe-i) (A; — 1)) G "Di estimates zi^-i by doing the intersection of 

sets d (Vi (k - 1)) and 5f } . Similarly to 03 Theorem 6] and using lfl4l Lemma 3], one can 
show that z^k-i = £i,fc-i with arbitrarily low probability of error if n is sufficiently large and 
R 1 <I{Y rl -Y l \U 1 ) + R\- ) . 

iii Using ^(ii fc-ilsi^-i) and Vi(k — 1), the receiver node X>i finally estimates the message i&i^-i = 
Wi if and only if there exists a unique codeword x™(w\) such that 

(2Ci(wi),«i(si, fc _ 1 ) l y5 1 (i - l),j/? l (ii J fc-i|si j fc_ 1 )) G One can show that u>i = iui, fc _i with 

arbitrarily low probability of error if n is sufficiently large and 

Ri < l(x 1 ;Y 1 ,Y r i I U,). (26) 



iv At the end of the block k, the relay looks for the suitable estimation of its observation that it 
intends to transmit to the receiver node V 1 by estimating z 1>k . It estimates z ljk = Z\ if there exists 
a sequence y^(zx\si t k) such that (w"(si ife ), y"i(^i|«si,fc), 2/" (&)) £ • There exists a such sequence 
if n is sufficiently large and R\ > I{Y r i\Y r \U\). 

From i, ii, iii we further obtain 

/(y^yi^y) < /(c^y). (27) 
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a 2 . [\mi\ 2 Re( Pi qZ) - {\ Pl \ 2 + Sl )Re(m in *)] + a r [K| 2 (|(ft| 2 + 1) - \rn\ 2 {\ Pi \ 2 + «<)] + {\q t \ 2 + l)Re{ mi n*) - \n t \ 2 Re( Pi q*) 



[\ Pl a r + qi\ 2 + Sid 2 , + l][\mia r + rii\ 2 + \ Pi a r + qt\ 2 + s.a 2 + 1] 



(28) 



The achievability proof for the second receiver node follows in a similar manner. Therefore, we have 
completed the proof. 
From the discrete case to the Gaussian case 

As mentioned in the beginning of this section, obtaining achievable transmission rates for Gaussian IRCs 
from those for discrete IRCs is an easy task. Indeed, the latter consists in using Gaussian codebooks 
everywhere and choosing the coding auxiliary variables properly namely choosing U\, U 2 , Y r>1 , and Y r;2 . 
The coding auxiliary variables U\ and U 2 are chosen to be independent and distributed as U\ ~ A/"(0, uiP r ) 
and U 2 ~ A/"(0, ^ 2 F r ). The corresponding codewords and u 2 convey the messages resulting from the 
compression of Y r . The auxiliary variables Y r i, Y r 2 write as Y r l = Y r + Zw} and Y r ^ 2 = Y r + Zwl where the 
compression noises Zwz ~ 

A/"(0, Nwz ) and Z^l ~ J\f(0, N^}) are independent. At last, the relay transmits 
the signal X r = U\ + U 2 as in the case of a broadcast channel except that, here, each destination also 
receives two direct signals from the source nodes. By making these choices of random variables we obtain 
the desired rate region. 

Appendix B 

Proof of Theorem \3A1 (optimal amplification gain for ZDSAF in IRCs) 

Using the notations given in Theorem 13 .41 and also the signal-to-noise plus interference ratio in the capacity 
function of Eq. (flOl) the rate Ri can be written as: 

i?j(fl r ) = C 

We observe that Ri(0) = C ^ |^p+r j an( ^ mat we nave an horizontal asymptote 

(1 771 ■ P \ 
- — ^ . Also the first derivative w.r.t. a r isa r is given in (1281) 
\Pi\ + S iJ 

The explicit solution, a* depends on the channel parameters and is given here below. We denote by 
A the discriminant of the nominator in the previous equation. If A < 0, then in function of the sign of 
\mi\ 2 Re(piq*) — (\pi\ 2 + Si)Re(min*), the function Ri(a r ) is either decreasing and a* = or increasing and 
a* = a r . Let us now focus on the case where A > 0. 





mid 




2 


\Pid r + qi 


2 + SiO 2 . + 1 
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1) If |m i | 2 Re(p i g*) - (|p;| 2 + s i )Re(m i O > then 

a) if c^l < and of} < then a* = a r ; 

b) if eg > and eg < then 

i) if a r > eg then a r * = 0; 

ii) if a r < eg then 

• if Ri(0) > Ri(a r ) then a* = else a* = a r ; 

c) if eg < and eg* > then the analysis is similar to the previous case and a* £ {0,a r } 

(2) 

depending on Or this time; 

d) if eg > and eg > 

i) if eg < eg 

A) if a r < eg then a* = a r ; 

B) if eg < a r < eg then a* = eg; 

C) if a r > eg then 

. if i?;(cg) > i?i(a r ) then a* = eg else a* = a r ; 

ii) if eg > eg then the analysis is similar to the previous case, exchanging the roles of eg 
and eg; 

iii) if eg = eg then a* = a r . 

2) If |mj| 2 Re(pig*) — (\pi\ 2 + Si)Re(mjn*) < then the analysis follows in the same lines and a* £ 
{0, a r , c r i , c r i }. 

Appendix C 

Proof of Theorem 03] (existence of an NE for the bi-level compression EF protocol) 

We want to prove that for each user Fq 9 ' is concave w.r.t. 9 9 . Consider w.l.o.g. the case of user 1. 
The general case of complex channel gains is considered. We analyze the second derivative of R q given 
in Eq. flT]). For the sake of clarity we denote by N[ q) = \h rl \ 2 lM P r {9) + N[ q \ T = \h rl \ 2 u^P r {q) and 
T 1 = \h 21 \ 2 6 yq, P 2 + N± q '. After some manipulations we obtain the following relation: — = Mi — M 2 
with M fc = gg*, k £ {1,2} where (for the sake of clarity we have denoted h\f by h ij ):NM 1 , NM 2 , DM U 
DM 2 are defined by (Til,®. 

We observe that the terms A k > 0, k £ {2, . . . , 7}. Also we can easily see from Eq. (|29l) that M 2 > 0, 
DM\ > 0. Thus if we prove that NM\ < the concavity of will be guaranteed. In this purpose we 
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NM 1 



NM 2 

DM 1 
DM 2 



\h 11 \'P 1 2 \h lr \- 1 -A 2 2 P 1 2 )\h 11 \ 2 e[ q) Px 

r a 5 



A 8 |fci 



r a 5 



A 7 (|feii| 2 Pi 2 |fei,-| 2 -A 2 2 Pi 2 )r 1 



, 9 A 7 A_8 IV o 

+ 2 iVa 5 * ~ J 



A 5 2 r 



1 + At 

1 ^ A 5 ' 

DM 2 , 



\h, I 2 T\ P, I 1 A 7 A 8 ri 



(29) 



A 2 
A 3 
A 4 
A 5 
A fi 
A 7 
A s 



2Ke(h u h* lr h* 21 h 2r ), 
\h u hl r \ , 

\h u \ 2 9[ q) P 1 + T 1 , 



N, 



(?) 



2/)(?) 



r(«) 



|/i 2j .| 2 * ) p 2 + iv r (9) + ivS 1 , 

Ael/inl^^Px + rxl^l^PPi, 



|/i 11 | z P 1 A 3 + |/n r | z PiA 4 - 2A 2 f^ ) P 2 - Ai^ ; PiP 2 



3(9) 



(30) 



plug the expressions of A 5 , A 6 , A 7 , A 8 into Eq. ([29]) and obtain that NM 1 = gjJgjL with NNM 1 , DNMi 
given in CD}, dSU). 

d 2 P (9) 

Therefore we obtain the desired result NM\ < and thus Mi > 0, which implies that . ,\ 2 < 0, 
whatever the channel parameters. 

Appendix D 

Proof of Theorem 03] (Number of Nash equilibria for ZDSAF) 

Before discussing these situations in detail, let us first observe that the two functions Fi(9j) are decreasing 
w.r.t. 9j and also FAO) = A , FAQ?) = where 9° = 

J Ca J J Cij 

In this section, we will investigate the NE of the game and also their asymptotical stability of each NE 
point. A sufficient and necessary condition that guarantees the asymptotic stability of a certain NE point is 
related to the relative slopes of the best-response functions and is given by 11371 11381 : 



January 12, 2011 



DRAFT 



29 



DNM l 



!r| 2 + 



NNM, = 2P?T if (^ ) )Vi| 4 W a + ^(»? ) ) 2 |^i| a |'»ii| a |fc 

(7V 1 (9) ) 2 |/ i lr| 2 + 2P 2 iV 1 (9) ^ ) |/l 21 f% r | 2 

( I h 2l | 2 P 2 + T + A^) (d { 2 q) P 2 # |/* 21 | 2 + 1^1 2 ^ P 2 + rt q) N[ q) 



(\h lr h 21 - huhrf ei q) P 2 + \h lr \ 2 N[ q) + |fcn| 2 #) 



(31) 



> 



- \h u h 2r - h 2l h lr \ 2 e[ q) p 2 N[ q) - \h 2l \ 2 \h lr \ 2 e[ q) p^ p 2 N[ q) - \h n \ 2 e[ q) pM q) N[ q) - 

\h 21 \ 2 \h 2r \ 2 (e^) 2 PiN[ q) - 2 \h 2l f 9 {q) P 2 Nl q) N[ q) - \h 2l f (^ } ) 2 p|iV r (9) - 

\h 21 h lr - h u h 2r \ 2 \h 21 \ 2 e iq) p 1 [ef) 2 pi - \k\*o®p x (n?) 2 - \h 2T \ 2 e^p 2 (n® 

(n^N^ - \h 21 \ 2 6 {q) P 2 N^T - N^ q) T N[ q) - \h 2r \ 2 9^ P 2 N[ q) T - 
\h 11 \ 2 \h 21 \ 2 9 {q) P 1 ei q) P 2 N^ 



< 



(32) 



dBRi dBR 2 



< 1 



(33) 



d9 2 d6 1 

in an open neighborhood of the NE point. We denote by V(#i, 6 2 ) an open neighborhood of (9i,9 2 ) G [0, l] 2 . 
1) If dx < and d 2 < 0, then the BR are constants BR^) = and thus the NE is unique (9f E , 9^ E ) = 

o for (e h e 2 ) g v(o,o) 



dBRi dBR 2 
d6» 2 d(9i 



(0, 0), for all cu > 0, Cji > 0. The condition (1331) is met since 
and thus the NE is asymptotically stable. 
2) If d\ < and d 2 > 0, then it can be checked that the NE is unique, for all Cu > 0, > 0: 9f E 
and 

' ^ , if d 2 < c 22 , 

C22 Z ZZ, 

1 , otherwise. 



qNE 



It can be checked that 



dBRi dBR 2 
d6 2 d^i 



for 2 ) G V(0f E , 9^ E ) and the NE is asymptotically stable. 



3) If d\ > and d 2 < 0, then, similarly to the previous item, we have a unique NE, for all cu > 0, 
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Cjl > 0: 0^ E 



and 



NE 



^ , if 

en 

1 , otherwise. 



di < en, 



Here as well we have 



dBRi dBR 2 
d8 2 d0i 



for (9 U 2 ) G V(0f E , 0^ E ) and the NE is asymptotically stable. 
4) If d\ > and d 2 > 0, we have to take into consideration the parameters cu > 0, Cji > 0. 

a) If -Fi(l) > 1 and F 2 (l) > 1, then we have di > c\ 2 + en and d 2 > c 2 \ + c 22 . In this case the 
BR are constants i.e., BRj(0 3 -) = 1 and thus the NE is unique (0f E ,0^ E ) = (1,1). We have 

= o for (0i, 2 ) e V(l, 1) and the NE is asymptotically stable. 

b) If .Fi(l) > 1 and F 2 (l) < 1, then we have di > c i2 + c n and d 2 < c 2X + c 22 . Here also the NE 
is unique and 8f E = 1 and 



NE 



<h ^ L , if d 2 > c 22 , 

«22 Z ZZ ' 



dBRi dBR 2 
d6» 2 d6»i 



, otherwise. 

for (0i, 2 ) e V(0f E ,0^ E ) and the NE is asymptotically 



Similarly, we have 
stable. 

c) If -Fi(l) < 1 and F 2 (l) > 1, then we have d\ < c i2 + en and d 2 > c 2i + c 22 . Here also the NE 
is unique and 0^ E = 1 and 

' ^ , if d 1 > cn, 

"--11 

, otherwise. 



NE 



dBRi dBR 2 
d6» 2 d0i 



for (0i, 2 ) e V(0f E ,0^ E ) and the NE is asymptotically 



Here as well we have 
stable. 

d) If -Fi(l) < 1 and F 2 (l) < 1, then we'll have di < c i2 + c n and d 2 < c 2i + c 22 . This case is the 
most demanding one and will be treated in detail separately. 
At this point an important observation is in order. The discussed scenarios, for which we have determined 
the unique NE, have a simple geometric interpretation. If the intersection point (0i,0 2 ) is such that either 
0* e R \ [0, 1] or 9* e R \ [0, 1] then the NE is unique and differs from this point ((9f E , 0^ E ) ^ (0*, 0*)). 
The case 4.(d) corresponds to the case where the intersection point (0j|\0 2 ) G [0, l] 2 is an NE point. Now 
we are interested in finding whether this intersection point is the unique NE or there are more than one NE. 
If < d± < cn + c i2 and < d 2 < c 22 + c 2 \ we have the following situations: 
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1) If C11C22 = C21C12, then the curves described by 0j = Fi(9j) are parallel. 

a) If d± — d-i, then the curves are superposed. In this special case we have an infinity of NE that 
can be characterized by (0f E , 0^) G T where: 

T={(9 1 ,9 2 )e[0,i} 2 \9 1 = F 1 (9^ E )}. 



dBRi dBR 2 
d6» 2 dfli 



1 for (0 l5 2 ) G V(0? E ,0^ 



In this case we have an infinity of NE such that 
and the NEs are not stable states. This can be easily understood since a small deviation from a 
certain NE drives the users to a new NE point. Thus, the users don't return to the initial state, 
b) If d\ 7^ di, then the two lines are only parallel. In this case it can be checked that the NE is 



unique and also asymptotically stable since again ^-^^ = for (0 1; 2 ) G V(0f E , 0* E ). In 
order to explicit the exact relation of the NE, one has to consider all scenarios in function of the 
sign of the following four relations Fj(0) — 1 and 9° — 1, % G {1, 2}. We will explicit only one 
of them. Let us assume that Fj(0) — 1 < and 9° < which means that di < min{ci 2 , £j ^f 1 } 
and d 2 < min{c 2 i, c 22 }. Here we have two sub-cases: 

. if A < A then the NE is characterized by 9f E = and 0? E = 

C12 C22 ^ 1 * C22 

. If A > A then the NE is characterized by 0^ E = ^23. and #ne = g 

C12 C22 J 1 C12C21 z 

2) Consider cuc 22 7^ C2i c i2- Here we have to consider all cases in function of the sign of the four 
relations Fj(0) — 1 and 9° — 1, i G {1,2}. We will focus on only one of them. Let us assume that 
Fj(0) — 1 < and 9° — 1 < and thus d x < min{ci 2 , c n } and d 2 < rnin{c 2 i, C22}. Here we have four 
sub-cases: 

. if A < A an d A > A then the NE is un j que: #ne = Q * md ^ne = Q * Also we haye that 

C 2 2 Cl2 Cn C 2 1 ' 1114 4 

< 1 for (0i, 2 ) G V(0*,0 2 ) and the NE is asymptotically stable. 



dBRi dBR 2 
d6» 2 d9i 



. If ^ > ^ and 

C22 C12 

//)NE ziNE^ 



'1 



— < — , then there are three different NE: 

Cll C21 

rf 2 



#2 ) e {(#*> #2); (0) (|S the intersection point and two other NE's on the border. 



The intersection point is unstable since 
two NE's are asymptotically stable since 
V(£,0). 
. If ^ 



dBRi dBR 2 
d6» 2 d9i 

dBRi dBR 2 
d0 2 d0i 



> 1 for 
= for 



2 ) G V{91,9*) but the other 



M2) eV(O^) and (0 l5 2 ) G 



and ^ < ^, then there are only two different NE: (0f E 0^ E ) G {(0, 



A 

C21 : 



A" 

C22 - 



di 



L OH 

C22 C12 " ' Cll C21 7 " " J - ' \- 1 J - 4 I L \ 1 C22 ' ' V Cn ' ' J ' 

In this case both of NEs are on the border, one of which represents the intersection point of 



the BR's. It turns out that the intersection point is not a stable NE because 



dBRi dBR 2 
d6» 2 d0i 



> 1 for 
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(0i,02) G However, the other NE is asymptotically stable since 

(0i,0 2 )eV(^,O). 

. if > A. anc i A. = A > t h en there are two NE: (0f E ,0P) G {(tS 0), (0, -^)}. Here the 

C22 ^-12 ^11 ^21 C22 

analysis of the stability of the two NE's is similar to the previous case. 
In conclusion, the number of NE states depends on the geometrical properties of the best-response 
functions. Three different cases can be identified: 1) when the lines 0j = Fi(9j) are superposed the game 
has an infinity of NE which are not stable; 2) when the lines have a unique intersection point that lies 
outside of the borders [0, 1] x [0, 1], the NE is unique and asymptotically stable; 3) when the lines have a 
unique intersection point (0^, 0|) that lies inside [0, 1] x [0, 1], there can be one, two or three different NE 
among which one is identical to this intersection point. In the case where the the NE is unique, it is also 
asymptotically stable. When the game has two or three NE, the intersection point (0^,0^) is an unstable 
equilibrium while the other/others are asymptotically stable. The best-response algorithm converges to one 
of the NE points depending on the initial state of the system. 
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Fig. 1. System model: a Q-band interference channel with a relay; q is the band index and q £ {1, Q}. 
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Fig. 3. Achievable system sum-rate versus x r (abscissa for the relay position) for a fixed y r (y r = 0.5do), with AF, DF and bi-level EF. 




Fig. 4. Best replies for a system composed of an IC in band (1) and IRC in band (2) when the ZDSAF protocol is assumed (fixed amplification 
factor). The number of equilibria is generally three as indicated the figure. 
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Fig. 5. ZDSAF relaying protocol with fixed amplification gain. Achievable network sum-rate at the NE as a function of A r £ [0, a r ] for 
L = 10m, e = 0.5m, Pi = 20dBm, P 2 = 23dBm, P r = 22dBm, Ni = lOdBm, N 2 = 9dBm, N r = 7dBm, 7 (1) = 7 (2) = 2. The optimal 
amplification gain A* = 0.05 < a r (l, 1) = 0.17 meaning that saturating the relay power constraint is suboptimal. 




Fig. 6. ZDSAF relaying protocol, full power regime. L = 10m, e = lm, Pi = 20dBm, P 2 = 17dBm, P r = 22dBm, Ni = lOdBm, 
N 2 = 9dBm, N r = 7dBm, 7 (1) = 2.5 and 7 (2) = 2. (a) Achievable network sum-rate at the NE as a function of (x-jz,y-R.) € \—L,L\ 2 (the 
optimal relay position (xtziVii) = ( — 1-2, 1-7) lies on the segment between Si and T>i). (b) Power allocation policies at the NE (8^ E , 8 2 E ) 
as a function of {x-r, yu) £ [— L, L] 2 (the regions where the uses allocate their power to IRC are almost non overlapping.). 
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Fig. 7. EF vs. DF relaying protocol. Achievable network sum-rate at the NE as a function of v £ [0, 1] for L = 10m, e = lm, Pi = 22dBm, 
P 2 = 17dBm, P r = 23dBm, TYi = 7dBm, N 2 = 9dBm, N r = OdBm, 7 (1) = 2.5 and 7 (2) = 2. The optimal relay PA v* = 1 is in favor 
of the better user and outperforms the uniform relay PA v = 0.5 for both EF and DF. 
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